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Abstract 
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to  the  multinomial  distributions  are  extended  to  more  distributions. 
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1 .  Introduction . 

Concepts  of  positive  dependence  of  sets  of  random  variables  (rv's) 
have  received  a  lot  of  attention  recently.  Their  study  was  found  to 
yield  a  better  understanding  of  the  structure  of  some  widely  used  multi¬ 
variate  distribution  functions  (df*s).  In  addition  to  this,  various 
useful  inequalities  were  obtained  with  applications  in  many  areas  of 
probaLility  and  statistics.  Barlow  and  Proschan  (1975),  Ch.  5,  include 
a  review  of  most  of  the  work  done  prior  to  1972.  A  list  of  more  recent 
references  can  be  found  in  Ahmed  et  al.  (1978). 

On  the  other  hand  notions  of  negative  dependence  have  received  very 
little  attention  in  the  literature.  Some  negative  dependence  analogs  of 
positive  dependence  concepts  have  been  mentioned  by  some  authors  (Lehmann 
(1966),  Brindley  and  Thompson  (1972),  Dykstra,  Hewett  and  Thompson  (1973) 
and  Shaked  (1977)  among  others).  In  the  bivariate  setting  the  random 
vector  ^j*^)  is  usually  said  to  satisfy  some  negative  dependence 
condition  if  (T^,-^)  satisfies  the  analogous  positive  dependence 
condition.  However,  this  method  of  formulation  cannot  apply  to  higher 
dimensions.  To  the  best  of  our  knowledge  Lehmann  (1966)  in  the  bivariate 
setting  and  Mallows  (1968)  in  the  multivariate  .setting  came  the  nearest 
to  a  systematic  study  of  negative  dependence  concepts;  Mallows'  discussion, 
however,  is  restricted  to  the  multinomial  distribution.  The  development 
below  is  in  the  spirit  of  Mallows;  his  results  are  special  cases  of  ours. 

While  the  first  draft  of  this  paper  was  being  written  two  related 
works  were  brought  to  our  attention.  The  first  work  by  Ebrahimi  and  Ghosh 
(1980)  discusses  some  negative  dependence  analogs  of  well  known  positive 
dependence  concepts.  Some  of  our  definitions  overlap  those  of  Ebrahimi  and 
Ghosh  (1980);  however,  our  main  results  differ  from  theirs.  The  second 
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related  paper  Is  by  Karlin  and  Rinott  (1980).  They  introduce  a  negative 
dependence  notion  which  is  closely  related  to  one  of  ours  and  they  obtain 
some  results  which  are  similar  to  ours.  Some  remarks  about  the  relationship 
between  these  two  works  and  the  present  paper  will  be  given  throughout  the 
text . 

The  main  motivation  for  our  definitions  is  to  try  to  formulate  the 
intuitive  requirement  that  if  a  set  of  negatively  dependent  random 
variables  is  split  into  two  subsets  in  some  manner  then  one  subset  will 
tend  to  be  'large'  when  the  other  subset  is  'small'  and  vice  versa.  In 
Section  2  we  define  the  conditions  to  be  discussed  We  derive  some  in¬ 
equalities  in  Section  3  and  prove  some  preservation  properties  in  Section  A. 
Examples  are  given  in  Section  5. 

In  the  following  "increasing"  stands  for  "nondecreasing"  and  "decreasing" 
for  "nonincreasing".  Vectors  in  lRn  are  denoted  by  t  •  (t, ,...,t  )  and 
t  <  means  t ^  <  tj,  i  *  l»...,n.  Similarly  £  <  t/  means  t^  <  tj, 

I  *  l,...,n,  and  0  «=  (0,  •••,())•  A  real  function  on  Hn  will  be  called 
increasing  if  it  is  increasing  in  each  variable  when  the  other  variables 
are  held  fixed. 

A  rv  X  is  said  to  be  stochastically  smaller  than  the  rv  Y  (denoted 

by  X  Y)  if  P(X  >  x)  <  P (Y  >  x)  for  every  real  x.  The  random  veccor 

X  *  (X, . X  )  ia  said  to  be  stochastically  smaller  than  Y  ■  (Y, , . . . ,Y  ) 

—  in  —  in 

[denoted  by  X  V  Y]  If  g(X)  V  g(Y)  for  every  g  c  C  where  C  is  the 

class  of  Borel  measurable  Increasing  functions  on  lRn.  If  X  and  Y  have 

St  St 

the  same  df  then  we  write  X  «  Y.  It  is  well  known  that  X  <  Y  if  and 
only  if 

P(XeU)  £  P(YeU)  for  every  upper  Borel  set  U  in  Tln. 


(1.1) 


[U  is  an  upper  set  if  x  e  U  and  x  <_  £  implies  that  £  e  U.]  According 
to  Kamae,  Krengel  and  O'Brien  (1977),  we  need  only  consider  open  upper 
sets  U  in  JRn.  It  is  also  well  known  that  for  every  random  vector  X^ 

(1.2)  X  +  m  V  X  whenever  a  () 
and  that 

(1.3)  X  +  A  ^  X  whenever  A  is  a  nonnegative  random  vector. 

Also  if  P(X  >  0)  -  1  then 

(1.4)  aX  V*  X  whenever  a  ^  1. 

See  Arias  and  Lehtonen  (1978)  for  an  excellent  review  on  stochastic 
ordering. 

If  X  and  Y  are  random  vectors  such  that  X  given  that  X  *  X. 
is  stochastically  smaller  (larger)  than  X  given  that  Y  ■  ,  whenever 
y  y*  and  £  and  y/  are  in  thc  support  of  Y,  then  we  write 

(x I V  -  y]  V  (?5 

More  precisely,  this  means  that  for  every  upper  set  U,  there  exists  a 
version  of  P €  U|Y)  -  <>(Y)  such  that  $(y)  is  increasing  (decreasing) 
in  y  on  the  support  of  Y. 

A  bivariate  function  K(*,*)  which  is  defined  on  x  (where 

and  are  subsets  of  JR)  is  said  to  be  totally  positive  of  order  2 

(TP2)  on  S  x  s2  if  K(x,y)  >_  0  and  if 

K(x,y)  K(x,,yt)  ^K(x,y’)  K(x’,y)  whenever  x  £  x* ,  y  ±  y* 


(1.5) 
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(see  Karlin  (1968)),  The  function  K  is  said  to  be  reverse  regular 
of  order  2  (RR2)  on  Sl  x  S2  K(x»y)  L  0  and 

(1.6)  K(x,y)  K(x'#yf)  <_K(x,yf)  K(x',y)  whenever  x£x',  y  <^y', 
(see  Karlin  (1968),  p.  12). 


2 .  Negative  dependence  concepts. 

Most,  but  not  all  of  the  positive  dependence  concepts  discussed  in  the 
literature  have  negative  dependence  analogs  that  can  be  obtained  by  changing 
the  direction  of  the  monotonicity  or  of  the  inequalities  which  define  them. 
Here  we  define  four  conditions  of  this  type  which  we  find  useful  because  we 
have  methods  of  identifying  distributions  which  satisfy  them.  Two  of  them 
are  direct  analogs  of  ’’conditionally  increasing  in  sequence”  and  ’’positive 
orthant  dependent”.  The  third  has  a  positive  dependent  counterpart,  but  it 
does  not  seem  to  have  been  discussed  in  the  literature,  while  the  fourth  is 
a  variation  of  "totally  positive  of  order  two  (TP^)  in  pairs”. 

In  the  case  of  positive  dependence,  one  of  the  strongest  and  most  use¬ 
ful  notions  is  that  of  TP^-ness  in  pairs;  that  is,  the  joint  density  or  the 
discrete  probability  function  f  is  assumed  to  exist  and  be  TP2  in  pairs 
(see  Barlow  and  Proschan  (1975),  p.  149).  The  natural  negative  dependence 
analog  then  is  to  assume  that  f  is  RR2  in  pairs.  There  are  several 
drawbacks  to  this  notion,  however.  Firstly,  unlike  the  situation  in  which 
the  joint  density  (or  discrete  probability  function)  is  TP^  in  pairs,  the 
marginal  densities  do  not  necessarily  enjoy  the  same  property.  A  simple 
3x2x2  discrete  example  suffices  to  show  this.  In  fact,  Theorem  5.1, 
p.  123  of  Karlin  (1968)  la  false  when  TP^  is  replaced  by  RR2.  (It 


5 


should  also  he  remarked  here  that  even  in  the  TP^  case,  one  must  make 
some  assumptions  on  the  nature  of  the  set  {f  >  0}  before  one  can  use 
the  above  Theorem  5.1  to  conclude  that  the  marginal  densities  are  also 
TP2  in  pairs:  see  Kemperman  (1977)).  A  possible  alternative  then  is  to 
assume  that  not  only  f,  but  all  its  marginal  densities  are  RR^  in  pairs. 
But  even  under  this  assumption  we  have  not  been  able  to  show  that  our 
weakest  condition  (see  Definition  2.3)  is  consequently  satisfied.  Ebrahlml 
and  fihosh  (1980)  claim  to  have  proven  this  result;  however ,  their  proof  is 
based  on  an  implication  which,  as  will  be  shown  below,  does  not  hold  (see 
discussion  after  Definition  2.4). 

Because  of  these  drawbacks  and  since  we  do  not  always  want  to  assume 

the  existence  of  a  density,  we  prefer  to  work  directly  with  the  measure 

itself.  This  point  of  view  is  consistent  with  our  other  definitions  in 

that  we  do  not  make  any  assumptions  on  the  existence  of  a  density  in 

defining  them.  We  are  thus  led  to  our  first  definition. 

Let  p  be  a  probability  measure  on  the  Borel  sets  in  JRn.  If 

I I  are  intervals  in  1r\  we  define  the  set  function  p(I~,.,.,I  ) 
in  in 

by  u(I-,...,I  )  ■  p(I  x...*I  ).  By  abuse  of  notation,  we  write  p  Instead 
In  In 

of  p.  If  I  and  J  are  intervals  in  ]r\  we  write  I  <  J  if  x  t  I, 
y  e  J  Implies  x  <  y,  that  is,  I  lies  to  the  left  of  J. 

2 

Definition  2.1.  Let  p  be  a  probability  measure  on  F  .  We  say  that  p 
is  reverse  regular  of  order  two  (RR2)  if 

(2.1)  M(I1,I2)  ji(IJ.IJ)  1  pdj.Ij) 

for  all  intervals  1^  <  Ij,  I2  <  I2  in  F1 .  We  also  say  that  p(IltI2) 
is  RR2  in  the  variables  I^,!^  **  1*  is  a  probability  measure  on  Fn 


(n  >_  2),  we  say  that  p  is  RR^  in  pairs  if  p(I^,...,In)  is  RR^  in 

the  pairs  1^,1  for  all  1  i  <  ]  <_  n  when  the  remaining  variables  are 

held  fixed.  The  random  variables  T.  t . . . fT  (or  the  random  vector  T  or 

in 

its  distribution  function  F)  are  said  to  be  RR^  In  pairs  if  its 
corresponding  probability  measure  on  TRn  is. 

Remarks: 

(i)  The  obvious  TP2  definitions  for  p  are  obtained  by  reversing 
the  inequality  in  (2.1). 

(11)  Clearly,  if  p  is  RR2(TP2>  in  pairs,  then  so  are  all  marginals. 
Furthermore,  it  is  not  difficult  to  show  that  if  F  is  the 
distribution  function  associated  with  p  and  if 
F(ti,...,tn)  -  p((t1,*>),  (t2, ,(tn,®))  is  the  survival 
function,  then  the  functions  F  and  F  are  RR2(TP2)  in  pairs 
in  the  sense  of  (1.5)  and  (1.6). 

(ill)  It  is  easy  to  show  by  a  simple  limiting  argument  that  if  p  is 
RR2(TP2)  in  pairs,  and  if  p  has  a  density  f  with  respect  to 
a  product  measure  m  *  m^*...xm^  oi  o-finite  measures  such  that 
f  is  continuous  on  the  support  of  m  and  zero  off  the  support 
of  m,  then  f  is  RR2(TP2>  in  pairs, 

(iv)  In  the  n  «  2  case,  we  have  the  stronger  converse,  namely,  if 

p  has  a  density  f  with  respect  to  a  product  measure  m  -  *  m 

of  a-finite  measures  which  is  RR2(TP2)  on  *  S2>  where 
is  the  support  of  tn^(i  *  1,2),  then  p  is  RRjfH^* 

(v)  In  the  TP2  case,  one  can  generalize  to  higher  dimensions  if  one 
makes  some  assumption  on  the  set  (f  >  0}.  Let  p  have  a 
density  f  with  respect  to  a  product  measure  m  *  TOi***%xmn  of 
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o-finite  measures.  Let  be  the  support  of  m^ .  Then 

the  support  of  m  is  S  =  S ,  x  . . . x S  .  We  assume  that  there 

In 

exists  S  =  S/..,xS  such  that  {f>0}f\S  =  S  and  that  f 

is  TP^  in  pairs  on  S.  Then  p  is  TP^  in  pairs.  Just 

use  Theorem  5.1,  page  123  of  Karlin  (1968)  repeatedly  to  show, 

e.g.,  that  for  fixed  intervals  I0,...,I  in  1R1, 

3  n 

g(x1,x2)  «  /.../  f(x1,x2>x3>..*fxn)  dm^Cx^) . . .dm^fx^) 

I*  I 
3  n 

is  TP 2  in  x^  and  x2  on  x  §2#  The  result  then  follows 
by  a  simple  integration  using  the  TP^ *  inequality  for  g. 

(vi)  The  generalization  to  the  RR2  case  is  not  as  simple.  If  one 
assumes,  however,  that  p  has  a  density  f  with  respect  to  a 
product  measure  m  *  m^x...*!^  of  o-finite  measures  such  that 
the  density  f  when  integrated  over  any  n-2  intervals  in  IR1 
is  RR2  in  the  remaining  unintegrated  variables,  then  p  is 
RR^  in  pairs.  In  terms  of  random  variables,  this  can  be  para¬ 
phrased  as  follows.  Let  T3,...,Tn  be  random  variables  with  a 
density  f  (with  respect  to  a  product  measure  of  o-finite 
measures) .  Then  p  is  RR2  in  pairs  if  and  only  if  for  every 
1  <  i  <  j  <  n  the  conditional  density  of 


(T.,T  )|  r\  {T.  C  I  } 

1  J  kjti.j 

is  RR2  in  and  tj  for  all  choices  of  intervals  I^(k^i,j) 

in  1R*.  Equivalently,  if  denotes  the  indicator  function  of 

I ,  then  p  is  RR^  in  pairs  if  and  only  if 


/.../  (  n 

1R  JR  k*i,j 


f  ( 


is  RR^  in  the  unintegrated  variables  t^  and  t^  for  all 

choices  of  intervals  I.  (k^i^j)  in  1R*.  By  replacing  XT 

Tk 

hy  4*^  in  the  above  integral,  and  requiring  it  to  be  RR^  in 
t^  and  tj  whenever  is  a  set  of  PF^  functions, 

one  obtains  the  negative  dependence  condition  of  Karlin  and 
Rinott  (1980),  It  is  a  condition  which  is  stronger  than  the 
one  of  Definition  2,1  as  can  be  easily  seen  by  recalling  that 
the  indicator  function  of  an  interval  is  PF2. 

(vii)  Clearly,  if  (y^)  is  a  sequence  of  RR^  in  pairs  probability 
measures  and  if  y^  converges  weakly  to  y  then  y  is  RR2 
in  pairs. 


Definition  2.2.  The  rvfs  T^ , . . , (or  the  random  vector  T  or  its 
df)  are  said  to  be  conditionally  decreasing  in  sequence  (CDS)  if,  for 
i  =  1,2, .. . ,n-l , 


(2.2) 


^Ti+llTl  "  **Tj  "  f  (t^,...,t^). 


Definition  2.3.  The  rv's  T,  ,...,T  (or  the  random  vector  T  or  its 
-  1  n  — 

df)  are  said  to  be  negatively  upper  orthant  dependent  (NUOD)  if  for  every 


(2.3) 


P(T  >  O  <_  1i1p(Ti  >  t1) , 


They  are  said  to  be  negatively  lower  orthant  dependent  (NLOD)  if  for  every  t^, 


(2 . 3 1 ) 


P (T  <  t)  <  7T  P(T,  <  t.). 
- -  i-1  i  -  i 
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When  n  *  2,  (2.3)  and  (2.3* )  are  equivalent,  but  not  when  n  3 
(see,  for  example,  Ebrahimi  and  Ghosh  (1980)). 

The  next  concept  has  a  natural  positive  dependence  analog,  however, 
we  are  not  aware  of  any  place  in  the  literature  in  which  it  has  been 
discussed . 


Definition  2.4.  The  rv's  (°r  t^e  random  vector  or  its 

df)  are  said  to  be  negatively  dependent  in  sequence  (NDS)  if,  for 
i  =  2,3, ... ,n, 

(2. A)  [(Tj . T±_i ) | Tt  =  t±]  Bf  ti. 

Often,  to  verify  (2. A),  one  can  find  it  easier  to  verify  that  for 
i  =  1, . . . ,n, 


(2. A') 


I  (Tx . Tt 


l’Ti+l* 


•V'Ti 


'i' 


Then  clearly  (2, A)  holds. 

We  now  investigate  some  of  the  relationships  among  the  various 
definitions.  First  note  that  NDS  implies  both  NUOD  and  NLOD  and  these 
implications  are  sharp.  To  see  this,  use  methods  similar  to  Barlow  and 
Proschan  (1975)  to  show  that  (2.4)  implies  for  i  *  l,...,n-l. 


(2.5) 


[  (T 


V' 


••Vi 


i+l 


ti+1l  ?  [(Tj_ . T^lT 


i+l 


whenever  tj^  <_  tj+^,  (Although  Barlow  and  Proschan  assumed  the  existence 
of  a  density,  a  modification  of  their  proof  works.)  But  from  (2.5)  it 


follows  that 
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P(T, 


C1 . Tn 


Cn> 


P  (Tj  > 


»T 


n-l 


t  1>  P(Tn 

n- 1  n 


Cn> 


i  p(T1  > 


*  ,Tn-2  >  Cn-2> 


;  p(t  >  t±) 

i-n-1 


i  ...  i  ”  P(T.  >  t  ) 
1-1  1  1 


which  proves  (2.3).  The  proof  of  (2.4)  =>  (2.3')  Is  similar. 

It  is  not  difficult  to  construct  an  example  which  shows  that 
CDS  NUOD  and  that  CDS  ^  NLOD  and  hence  CDS*#  NDS.  For  example,  let 
Pdj-l.Tj-l)  -  P(T1-l,T2-2)  -  P(T1-2,T2-2)  «  .1  and  P(T1-2,T2-1)  -  .7 
and  let  given  T^  -  1,  T2  -  1  be  degenerate  at  11,  given 

Tj  =  1,  T2  -  2  be  degenerate  at  1,  given  -  2,  T2  “  1  be 

degenerate  at  10  and  given  -  2,  T2  -  2  be  degenerate  at  1. 

Then,  clearly,  (T^Tj.Tj)  is  CDS  but  P(T  >  1,  Tj  >  1)  -  .7  >  .64 
-  P(Tj  >  1)  P(T3  >  1),  thus  (T  T2,T3)  is  neither  NUOD  nor  NLOD. 

Fbrahimi  and  Ghosh  (1980)  claim  that  CDS#NU0D;  the  example  shows  that 
this  is  not  the  case. 

Next  we  show  that  NDS  ^  CDS.  Let  (T^.T^T^)  take  on  values  on 
the  eight  vertices  of  the  unit  cube  such  that 

P(T]=1,  T2-0,  T3=l)  -  P(Tj-l,  T2-l,  T3-0)  -  .2  and  the  other  six  pro¬ 
babilities  are  .1.  It  is  easy  to  see  that  T^  and  T 2  are  independent 
and  that  [(T^T^lTj  -  0]  V  [(T^T^I^  -  1].  Thus,  (T^Tj.Tj)  is 
NDS .  But  P(T3  >  0 1 Tx  -  0,  T2  -  0)  -  1/2  <  2/3  -  P(T3  >  -  1,  Tj  -  0), 

hence  (T^ *^2*^3^  CDS. 

We  will  shortly  show  that  under  some  reasonable  assumption,  RR^^CDS. 
It  is  not  known  whether  NDS,  but  it  does  imply  NUOD  and  NLOD.  This 

is  clear  since  from  the  remark  (ii)  following  Definition  2.1,  we  have  that 


.. nvTr  T'T* • 
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both  K  and  F  are  RR^  in  pairs.  The  result  then  follows  by  a 
simple  argument. 

The  implications  RR2=*CDS,  RR2=?NUOD,  RR^  NLOD,  NDS^NUOD 
and  NDS  NLOD  justify  the  consideration  of  the  RR2  and  the  NDS 

concepts.  In  Sections  4  and  5  it  will  be  shown  that  many  df*s  are  RR2 
or  NDS  and  thus  these  df's  satisfy  the  meaningful  CDS  concept  and  the 
inequalities  which  can  be  derived  from  the  NUOD  and  NLOD  concepts.  In 
addition  to  it  the  NDS  concept  is  intuitively  meaningful  by  itself.  We 
mention,  in  passing,  that  since  Karlin  and  Rinott's  (1980)  condition 
implies  the  RR2  condition  it  follows  that  it  implies  the  NUOD  and  the 
NLOD  inequalities  (2.3)  and  (2.3*).  Actually,  Karlin  and  Rinott  (1980) 
obtained  some  additional  useful  inequalities  which  follow  from  their 
stronger  condition. 

To  justify  calling  (2,l)-(2.5)  "conditions  for  negative  dependence" 
we  have  to  show  that  they  imply 


(2.6) 


C0V(TifTj)  <^0,  i  <  i  <  ]  <  n, 


when  the  second  moments  exist. 


From  (2.3),  it  follows  that  P(Ti  >  >  t^ )  <bP(T1  >  t±)  P(Tj  >  tj) 

and  it  is  well  known  that  this  inequality  implies  (2.6)  [see  e.g.,  Lehmann 
(1966)].  Similarly  (2.3*)  implies  (2.6)  and  hence  also  (2.4)  and  (2.5) 
imply  (2.6). 

If  (T., . . .  ,T  )  Is  CDS  then  [T2 1 T.  -  tj]  hence  covO^,^)  <  0. 


Thus,  if  (Tff(ij  , . . .  *T7r(n)^  is  CDS  for  every  permutation  tt  of  (1,2,.  ..,n} 
then  (2.6)  holds. 

We  close  this  section  with  a  proof  that  RR2  in  pairs  implies  CDS  under 
some  reasonable  assumptions.  Let  p  be  RR2  in  pairs  and  let  (T^,...,T^) 
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f >t*  any  random  vector  having  y  as  its  induced  probability  measure. 

Kix  t  l  IK  and  an  integer  i,  1  <_  i  <_  n-1  .  Define  the  set  functions 
v  and  A  on  the  Borel  subsets  of  TR*  by 

v  (A)  **  P  ^  ^14*1  *  * ' '  *  ^  ^ 


A (A)  -  P((T1,...,T1)  e  A). 

X  oo 

Note  that  A  is  just  a  marginal  of  y.  Then  if  is  a  sequence 

of  rectangles  which  partition  IR*  for  each  l  and  whose  mesh  size  tends 
to  zero  as  l  -*■  °°,  and  if  (R^  }  is  a  refinement  of  for  aH  * 

then  it  is  well  known  by  martingale  arguments  (see,  e.g.,  Meyer  (1966)) 
that  on  a  set  I)  with  A(DC)  *  0,  v^ik^*^Rik^  **  4>  (t  ^ ,  •  .  .  ,  t^,)  as 
decreases  to  (t^,...,t^)  pointwise  and  in  L1 ;  moreover  <f>  (t^  , .  •  • , 

is  a  version  of  p(Ti+1  >  | T1  =  . *  t^).  We  may  assume  without 

loss  of  generality  that  D  c  support  A.  We  want  to  show  that  <J>  is 
decreasing  on  D.  This  will  follow  if  we  can  assume  that  the  support 
A  of  A  satisfies  a  chain  condition;  that  is,  if  t  e  A  and  t*  e  A  with 


t  <  tf,  then  there  exist  t,  e  A  such  that  t  *  <  t-  < • • •<  t  “  tf  and 

— -  — j  —  — u  —x  - m  — 

t,  differs  from  t..,  in  only  one  component.  To  see  this,  suppose  that 

”1  “i+i 

t^  and  differ  in  only  the  first  component.  Then  let 

'i  *  !i+i  ■  ‘--■•w '  'i+i  ■  ‘w’  "d  h . ’i  be  any 

intervals  such  that  t^  e  1^  *  e  1^  x  If  P 

is  the  marginal  of  (Tj . T  j) ,  then  since  4t  is  RR2  ln  Palrs*  we  have 
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or  that 

^  (1  j  *^2  9  *  *  * 

A(i1,i2,...ti7)  -  H^T^TTTTTT^)  ' 

Note  that  there  is  no  trouble  dividing  since  the  denominators  are  non  zero. 
Iterating  this  procedure  to  pass  from  t.  to  t0,  to  .  to 

t  and  then  letting  the  intervals  shrink  we  find  that  $(t)  i  “Kt').  Now 


4>*(t) 


4>(t)  if  t  e  D 


inf  {<J>  (s) ;  £  e  D,  £  <_  t_}  (inf  0  -  +°°)  if  t  ^  D, 


It  easily  follows  that  4>*  is  decreasing  everywhere  and  is  a  version  of 

P(Ti+l  >  ci+llTl  “  h’*"’1!  =  ci)* 

The  chain  condition  is  easily  seen  to  be  satisfied  if  the  support 
of  vj  is  a  cross  product.  Without  some  type  of  chain  condition,  we  can 
only  show  that  $  is  decreasing  componentwise  on  D.  In  this  case,  one 
may  not  be  able  to  extend  $  such  that  it  is  decreasing  everywhere  and  is 
still  a  version  of  the  condition  probability. 


3.  Some  inequalities. 

This  section  is  devoted  to  the  derivation  of  some  inequalities  which 
may  be  of  special  interest.  The  results  closely  parallel  those  found  in 
Karlin  and  Rinott  (1980),  but  are  derived  under  weaker  assumptions. 


Proposition  3.1.  )  is  NUOD  if  and  only  if 


E[  "  MT<>]  i  "  Eft.  <T  )| 

i-1  11  i-1 


(3.1) 


J* 
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whenever  all  <\>^  are  nonnegative  and  increasing.  The  result  (3.1) 
is  also  true  if  we  replace  NUOD  by  NLOD  and  increasing  by  decreasing. 


Proof.  Let  4>.  (t)  ■  X.  v(t),  i  *  l,...,n.  Then  (3.1)  reduces  to 
the  NUOD  inequality  (2.3).  Since  each  side  of  (3.1)  is  multilinear  in 
the  the  result  holds  for  nonnegative  linear  combinations  of  such 

indicator  functions  and  hence  for  the  general  4^  by  a  standard  limiting 
argument.  |  j 


Now  suppose  that  p  is  RR^  in  pairs.  For  every  i,  let 
all  intervals,  with  <  K^. 


Theorem  3.1.  If  1  <_  k  <  n ,  then 


<3-2>  »<J1 . V  »«1 . V  i“<Ji . Wi . V 


*  U<I1 . Jk,Jk+l . Jn'' 


** 


The  result  (3.2)  Is  also  true  if  we  replace  all  J's  by  K's. 


Proof .  We  proceed  by  induction.  If  n  *  2,  then 


r 


w(JlfJ2)  ll(IltI2)  1 


w(lltJ2)  if  k-1. 

if  k  -  2. 


The  case  k  ■  1  follows  from  the  RR2  assumption  since 


li(J1,J2)  w(JltK2) 

0  > 

- 

p 

u(K1,J2)  wO^.Kj) 
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and  the  case  k  =  2  is  an  identity* 

Now  suppose  that  (3*2)  is  true  whenever  v  is  a  probability  measure 

on  lRn  which  is  RR^  in  pairs.  Let  u  be  a  probability  measure  bn  IRn+^ 

which  is  RR^  in  pairs  and  let  1  <  k  <_  n+1 .  Since  there  is  nothing  to 

prove  if  k  =  n  +  1 ,  we  may  assume  that  1  <_  k  <_  n.  Similary,  we  may 

assume  that  ,J  . , )  4  0.  It  suffices  then  to  prove  that 

i  n  n-ri 


(3  3)  m(I1 . In,In+l) _  <  ti(Il,,,,,In,Jn+l) _ 

u(J1,...,Jk,ik+1,...,in.in+1)  -  v(J1,...,Jk,ifcfl,...,in,Jn+1) 

since  by  the  induction  hypothesis f  we  have  that 

u(Ilt...,In,Jn+1) _  p(Ilt...,Ik>Jk+1 . Jn,Jn+l) 

p(J1,...,Jk,Ik+1 . In'Jn+l)  ~ 

But, 

u(Il,,,,,In,In-fl)  ti(J1,I2,...tIn,In+1) 
u(Il . In,Jn+l)  ”  m(J1,I2 . In,Jn+l^ 

p(J1,...,Jk,Ik+1,...>In,In+1) 
p(Jl . Jk  *  Xk+1  *****  Xn  *  Jn+l7 

which  is  another  way  of  writing  (3.3).  The  Jth  inequality  above  follows 
from  the  fact  that  y  is  RRj  in  the  pair  j  and  (n+1) .  |  | 

Suppose  that  (T-,...,Tn)  is  RR2  in  pairs. 

Corollary  1.  If  a,  $  partition  {l,...,n)t  then 
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(3.4)  P(T1eJ1,iea^8)  P(T  el^ieavr  8) 

<  P(TieJ1,lGo;TJcIj,je8)  P^elj.ieajTjCJj  ,JeB). 

It  also  holds  true  If  we  replace  all  J's  with  the  K's. 


Corollary  2.  If  a,  8,  Y  partition  {l,...,n},  then 

(3.5)  P(TieLi,ieot;TjeJj,j£6uY)  y) 

1  P(TieLi,iea;TjeJj  , j eB ^cl^key)  P(TieLi>iea;TjeIj  . jcB;TkeJk#kcy) 

for  any  intervals  L^,  i  e  a.  It  also  holds  true  if  we  replace  all  J's 
with  the  K's. 

If  we  take  »  (“ot>,bi]  and  1^  «  (-«>,“>)  for  i  €  av  B  in 
Corollary  1,  we  get  (3.5)  of  Karlin  and  Rinott  (1980).  If  we  take,  in 
Corollary  2,  *  [a^fb^],  i  e  a,  -  (-°°*bj  ]  and  1^  -  (-*0,®)  for 

1  e  Buy  we  get  (1.7)  of  Karlin  and  Rinott  (1980). 

Also,  note  that  as  soon  as  we  have  an  inequality  of  the  form 


<3-b>  P<T1  i  »1 . T„  i  b„>  i  p<Tl  i  bl . Tk  i  V  p<Tk«  i  bk+l . T„  S.  V 

it  follows  as  in  the  proof  of  Proposition  3.1  that 


(3.7) 


E[  *  4> .  (T. )]  <  E[  it  *.(T  )]  E[  tt  *.  (T  )] 

i-1  1  1  i-1  1  1  j-k+1  i  i 


whenever  $  are  nonnegative  and  decreasing.  Similarly,  if  we  have 


P(T1  5  al . Tn  >  an)  -P(Tl  51  al . Tk  >  V  P(Tk+l  >  ak+l,,,,,Tn  *  ®n 


(3.8) 
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then  (3.7)  holds  for  nonnegative  and  increasing. 

In  particular,  (3.6)  holds  if  F  is  RR^  in  pairs  and  (3.8)  holds 
if  F  is  RR2  in  pairs.  In  fact,  if  F(or  F)  is  RR2  in  pairs,  then 
it  is  MRR2  in  the  sense  of  Karlin  and  Rinott  (1980);  i.e., 

F(*  AZ.)  1F(x)  F(^). 


4 .  Closure  Results. 

Preservation  theorems  are  useful  for  identifying  negatively  dependent 
df’s  or  for  constructing  new  negatively  dependent  df's  from  known  ones. 

In  this  section  we  discuss  some  preservation  results  and  describe  a  method 
for  the  construction  of  negatively  dependent  df*s. 

Theorem  4.1.  If  T, ,  • . » ,T  are  (*)  and  if  •  are  strictly 

—  in  in 

increasing  functions  then  ^(T^) , . . .  •♦n0'n)  are  (*)  where  (*)  is 
one  of  the  following:  RR2  in  pairs,  CDS,  NDS ,  NUOD  or  NLOD. 

Theorem  4.2.  If  (T^#.,*fT  )  and  (S^,...,Sn)  are  independent  and  are  (*) 
then  (T^ , . . . ,T^ ,  is  (*)  where  (*)  is  the  same  as  in  Theorem  4.1 

The  proofs  of  these  theorems  are  straightforward  and  will  be  omitted. 
The  following  preliminaries  are  needed  for  the  statement  of  Theorem  4.3 

A  univariate  density  f  is  said  to  be  a  Polya  frequency  function  of 
order  2  (PF2)  f(x-y)  is  TP2  on  H  *  J*.  A  probability  function  f 

is  PF2  if  f(x-y)  is  TP2  on  JN  *  JN  where  K  -  {  •  •  •  A 


thorough  discussion  of  PF2  densities  and  many  examples  can  be  found  in 
Karlin  (1968). 
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It  will  ho  shown  in  a  later  section  that  many  multivariate  random 
vurinl  1  os  satisfy  a  certain  structural  condition.  In  the  next  theorem 
wo  will  show  that  this  condition  implies  some  of  the  dependence  conditions 
that  wo  have  introduced. 

Theorem  4.3.  Let  Sq,S^,...,S^  be  independent  rv's  and  assume  that 

each  has  a  PF^  density  (or  probability  function).  Fix  s  and  let 

(T- , . .  .  ,T  )  have  the  same  joint  df  as  the  conditional  df  of  ) 

L  n  in 

given  that  +  S-+...+S  *  s,  that  is, 

()  i  n 

(4.1)  (T1,...,Tn)  [(S1,...,Sn)|S0  +  S1+...+Sn  -  s]. 

Then  (T^,...,Tn)  is  RI^  in  pairs  and  consequently  CDS,  NUOD  and  NLOD. 

Proof.  Let  \i  be  the  probability  measure  of  (T.  , . . .  ,T  )  on  lRn.  Then 
-  1  n 

by  assumption  p(I_,...,I  )  *  P(S.  e  I-,...,S  e  I  I  S_  +  S.+...+S  ■  s) . 

1  n  ii  nnUi  n 

Now  the  loint  density  of  )  I  +  S-+...+S  “  s]  is  given  by 

i  n  u  i  n 

n 

c  tt  f  (s  )f  (s  -  s.-.  ..-8  ) 

1*1  1  1  0  1  n 

where  c  is  a  normalizing  constant.  We  first  show  that  p  is  RR2  in 
the  variables  I ^ ,  1^  when  the  remaining  intervals  are  held 

fixed.  According  to  the  remark  (vi)  following  Definition  2.1,  we  need 
only  show  that 

g(s1,s2)  -  c  f2(s2)  f0^S-8l",82^ 

is  in  s^  and  S2*  where 
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T 


111^*11  Ml'l  ■  •'■'I'  ,  ^ 

■Ti  1  .. 


■)  =  /•••/  f.(s  )...f  (s  )  f_(r-s  -...-s  )  dm(s  )...dm(s  ) 

<•  1  1  3  3  n  n  0  3  n  3  n 

I,  1 

)  n 

and  m  is  either  the  Lebesgue  measure  or  the  counting  measure.  However, 

the  above  is  nothing  but  the  convolution  of  the  PF^  functions 

XT  ,  ...,f  XT  and  fn,  where  X  is  the  indicator  function  of  the 
3  I  ~  n  I  U  A 

3  n 

set  A,  and  so  is  PF2*  It:  eas*-ly  follows  then  that  g  is  RI^ . 

The  proof  that  u  is  RR^  in  the  variables  1^,  1^  for  all  1  £  i  <  j  <_  n 
is  simi lar .  | | 

Remark.  Actually,  it  is  not  difficult  to  show  that  under  the  assumptions 

of  Theorem  A. 3  the  df  of  the  random  vector  (T  , . . .  ,T  )  satisfies  the 

1  n 

stronger  condition  of  Karlin  and  Rinott  (1980),  The  consequences  of  this 
oh? ervation  are  discussed  below. 

It  is  well  known  (see  Section  5)  that  the  multinomial,  the  multivariate 
hypergeometric,  the  Dirichlet  and  the  Dirichlet  compound  multinomial  random 
vectors  as  well  as  some  multivariate  normal  random  vectors  with  nonpositive 
correlations,  can  be  represented  as  in  (A.l).  Thus,  all  the  respective  df's 
are  RR^  in  pairs.  Karlin  and  Rinott  (1980)  have  shown  that  these  df’s 
actually  satisfy  their  stronger  condition;  however,  their  proofs  are  quite 
involved  and  differ  from  one  case  to  another.  By  the  Remark  after  Theorem  A. 3 
these  results  of  Karlin  and  Rinott  (1980)  follow  at  once. 

The  following  theorem  will  be  found  useful  in  the  next  section. 

Theorem  A. A.  Assume  that  (T.  , .  , .  ,T  )  and  (S, ,...,S  )  are  independent 
-  in  in 

and  NDS .  If  all  the  univariate  marginal  densities  (with  respect  to  Lebesgue 
measure)  or  probability  functions  in  the  discrete  case  of  and  T  are 
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! 
! 

? 

PF„ ,  then  (T,  +  S,,...,T  +  S  )  is  NDS. 

l  linn 

Remark .  Karlin  and  Rinott  (1980)  have  proven  a  similar  result.  They 
assumed  that  and  T  satisfy  their  RR^  condition  and  that  they 

have  PF^  marginals,  and  they  showed  that  then  £+  T  satisfy  some 
inequalities  that  are  essentially  variants  of  the  NUOD  and  the  NLOD 
inequalities . 

For  the  proof  of  Theorem  4.4  we  need  the  following  lemmas. 

Lemma  4.1.  Let  X  and  Y  be  independent  rv’s  and  assume  that  Y  has 

a  PF^  density  (probability  function).  Then 

[X | X  +  Y  =  z]  V  z. 

This  lemma  is  actually  Example  12  of  Lehmann  (1966). 

Lemma  4.2.  Let  X  and  Y  be  independent  and  assume  that  X  and  Y 

have  PF^  densities  (probability  functions).  Then 

(A. 2)  [ (X ,Y) | X  +  Y  -  z]  z. 

Proof.  Let  z1  <_  z^.  Denote  by  ^  the  upper  set  {(s,t):  s  >  x,  t  >  y}. 

First  it  will  be  shown  that  for  every  (x,y) 

(A. 3)  P{(X,Y)  e  U(Xfy)|x  +  Y  -  z^  <P((X,Y)  e  U(x,y)lX  +  Y  “  z2]  * 

Let  Aa  »  { (s,t) :  s  +  t  -  z± ,  x  <  s  <  z^  -  y} , 

B  ®  {  fs ,  t) :  s  +  t  -  z a  >zt  -  y),  Ct  -  {(s,t):  s  +  t  -  z^,  s  lx). 


mm 


i  =  1,2  (see  Figure  1).  From  Lemma  4.1  [x|x  +  Y  =  z]  St  z;  hence 

P{  (X,Y)  e  C1|X  +  Y  =  Zj }  =  P{X  <_  x|x  +  Y  =  z^ 

>  P{X  <  x|x  +  Y  -  z2>  -  P{  (X,Y)  e  C2|X  +  Y  -  z^ . 


Similarly 


P{  (X ,Y)  e  BjJx  +  Y  -  >_P{(X,Y)  e  B,,|x  +  Y  -  z.,) . 

Tli  us 

p{  (X ,Y)  e  Ax|x  +  Y  =  Zj}  >  P{  (X ,Y)  e  A2|x  +  Y  -  z^ 
which  is  (4.3). 

t 


Figure  1. 

m 

Now  consider  upper  sets  of  the  form  U  m  W  U ,  ,  called  funda- 

i-1 

mental  upper  domains  in  Block  and  Savlts  (1979) •  Without  loss  of  generality 
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assume  that  x  <  x,  <  ...  <  x  and  y1  >  yn  >  ...  >  y  .  Define 
1  ■  l  —  in  J  —  Z  —  —  m 


i ^  e  min  { i :  +  y^  z^} 

*  min  {i  >  i^s  +  y^  >  z 

and,  by  induction 

(A. 4)  iR+1  =  min  {i  >  jR:  +  y±  <_  z^) 


(4.5) 


Jk+1  “  mln  {1  >  ik+l: 


Xi  +  yi  >  Z2} 


k  =  l,2,...,n  where  n  <  00  is  the  largest  k  such  that  the  set  on  the 
right  hand  side  (RHS)  of  (A. 4)  or  (4.5)  is  not  empty.  If  there  are  n  i's 
hut  only  n-1  j's  define  j  *  m  +  1  (see  Figure  2  in  which  m  *  12, 

n  -  3,  i  »  3,  Jj  -  6,  ±2  -  8,  J2  -  11,  i3  -  12,  j3  -  13). 

-  n 


Let  U  -  \J  U 


(*i  ,yA  )’ 


In  the  following  the  first  inequality  follows 


k«l  9JA  J 

\  3k-l 

from  UC  U  and  the  second  one  from  (4.3): 


P{  (X ,Y)  e  U|x  +  Y  -  zL>  <  P{  (X ,Y)  e  U|X  +  Y  - 


n 

-  I  P{  (X ,Y)  €  u  Jx  +  Y  -  z  } 

k-1  ( W1 

(4.6) 

<  l  P{  (X,Y)  e  U  JX  +  Y  -  z  ) 

W1 


P{  (X ,Y)  e  U | X  +  Y  -  z2)  -  P{ (X ,Y)  e  U | X  +  Y  ■  z2). 
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Sin cv  open  upper  sets  can  be  approximated  by  fundamental  upper  domains  (Block 
anti  Savits  (1979))  it  follows  from  (4.6)  that  for  every  open  upper  set  V 

P{  (X,Y)  e  U | X  +  Y  -  Zj}  <  P{ (X,Y)  e  u|x  +  Y  -  *2) 

and  the  proof  of  the  lemma  Is  complete  by  (1.1).  || 

Lemma  4.3.  Let  X*  (X^,...,Xm)  and  Y_  =  (Yj , . . . .Y^)  be  independent 
and  assume 


(4.7) 


[  (Xj , . 


,X  .  )  |  X  -  x  ]  Y  X 
m-l  'mm  m 


and 


(4.8) 


[  (Yx , 


,Ym  ,)|y  -  y]  V  y  . 

m— i  mm  m 


Furthermore,  assume  that  X  and  Y  have  PF_  densities  (probability 

m  m  l 

f  unct ions) .  Then 


[(X,  +  Y  . X  +  Y  , )  I  x  +  Y  -  z]  ¥  Z 

I  i  m-l  m-l  m  m  m  \ 


Proof .  Clearly,  for  any  increasing  function  g, 

I.(p(X  +  Y X  +  Y  )|x  +  Y  «  z  ]  -  F[<KX  ,Y  )|X  +  Y  -  z  ] 

1  I  m-l  m-l  m  m  m  mm'm  m  m 

where  rp  (x  ,y  )  ■=  F[f»(X  +  Y  , . . .  ,X  ,  +  Y  ,)|x  -  x  ,  Y  -  y  ].  However, 

mm  l  i  m— i  m— I  m  mm  m 

-t>(x  ,y  )  decreases  in  x  and  in  y  because  of  (4.7),  (4.8)  and  in- 
m  m  m  m 

dependence.  Thus,  by  Lemma  4.2,  E  [  (X  ,Y  )|X  +  Y  ■  z  ]  decreases  in 

T  J  m  mimmm 

*  •  II 

m 


Proof  of  Theorem  4.4.  Let  1  e  {l,...,n-l}.  Substitute  m  ■  1  +  1  in 


2C> 


hemma  4.3  to  obtain 


l  +  S 


I’*’* 


V 


i+1 


Si+1 


=  h+l1 


V 


"i  +  1  ’ 


that  is,  T  +  S  is  NDS .  [  | 


*> .  Examples . 


r) .  1  .  Ttie  multinomial  df  . 

bet  )  have  the  joint  probability  function  with  parameters 

(N.pp  .  . .  .p^)  } 


P(T1  =  V 


,T 


tn)  * 


n: 


t  :...t  :(n-  l  t.): 

1  n  i 


n 

l 

i*l 
n 


n  ti 
<i*i pi  > 


N-  l  t 

n  .fr*  i  n 

*  (1-  l  Pi)  ,  t  >J,  l  t  <  N, 

i-1  1  i-1  1 


where  p  >  0  (i  =  l,...,n)  and  0  <  [  p  c  ], 

i-1  1 

The  multinomial  df  is  the  conditional  df  of  independent  Poisson 

rv’s  given  their  sum.  Thus,  by  Theorem  4.3  the  multinomial  df  is  RR^ 

in  pairs  and  hence  it  is  also  CDS,  NUOD  and  NLOD .  By  Remark  (iii)  the 

joint  probability  function  of  (T1f...,T  )  is  RR.  in  pairs.  By  the 

In  Z 

discussion  after  Theorem  4.3  the  multinomial  df  satisfies  the  RR^ 

condition  of  Karlin  and  Rinott  (1980). 

To  show  that  (T1lt,.,T  )  is  NDS  it  is  enough  to  show  that  for 
1  n 

0  <  t  <  N-l, 
n  —  * 


( (Tj 


r  .) |t 

n-l  1  n 


t+11 

n 


(5.1.1) 


T„-i>lT„  -  t„i  i  i<T, 


I'.M.mse  then,  by  symmetry  one  has  (2.4*).  The  left  hand  side  (LfIS)  of 

(3.1.1)  has  an  (n- 1 ) -d imensional  multinomial  df  with  parameters 

(N-i  n  ,<|  i  t .  .  .  #q  j  )  and  the  right  hand  side  (KIIS)  of  (5.1.1)  has  the 

same  df  with  parameters  (N-t  -1,  q,,...,q  ,)  where  q.  ■  p./(l-p  ). 

n  i  n-i  11  n 

Thus , 

(5.1.2)  LHS  (5.1.1)  S-  RHS  (5.1.1)  +  (S.,...,S  .) 

1  n- 1 

where  (Sj  ,  .  .  .  , S^_^)  has  a  multinomial  df  with  parameters  (lfqj,...q 

and  the  sum  on  the  RHS  of  (5.1.2)  is  of  independent  random  vectors. 

Inequality  (5.1.1)  follows  now  from  (5.1.2)  and  (1.3). 

het  (X.  )  he  th?  sum  of  independent  n~dimensional ,  not 

l  n 

necessarily  identically  distributed,  multinomial  random  vectors,  that  is 
let  (X^,...,X^)  have  the  probability  generating  function 

m  h 

(pUul+*  *  ,+PnfcUn) 

n 

where  pu  -  °*  a  “  l.*»«»n),  0  <  I  P1(l  <  1  and  >  1  (1  *1,2,. 

i*l 

The  univariate  marginal  df's  of  a  multinomial  df  Is  a  binomial  df . 

Since  the  binomial  df  has  a  PF^  probability  function  it  follows  from 

Theorem  4.4  that  (X. , . . . ,X  )  is  NDS.  Hence  (X-,...,X  )  is  also  NUOD 

In  in 

and  NLOI) ,  that  is, 

p<xi  >  . *„  >  *„>  iji  p<Ki  -  xi> 

and 

p<x,  —  xi . i «,»  i  til  p«,iV' 

These  inequalities  are  stated  unproven  in  Mallows  (1968).  Compare  also 
l.ehmann  (1966),  pp.  1  143,  1144  and  1151. 
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».2.  Multivariate  normal 

* 

Let  T  =  (Tp... 

•V 

he  i 

i  mu 

1 1  i variate 

svnmerric  normal 

rant!  rnn 

veetor  with  Corr(1\, 

V 

=  P  1 

o. 

1  i  v  <  1 

*  n .  Then  r>  > 

-  (n-l) 

Ive  will  show  that  T 

is 

rr2 

in 

pairs . 

Using  Theorem  4.1  assume*  without  loss  of  generality,  that  FT  ^  *  0 

and  Vnr(T  )  *  1,  i  =  l,...,n.  Let  Ylf...fY  he  independent  identically 
I  in 

distributed  normal  rv’s  such  that  KY^  =  0  and  Var(Y^)  *  1-p  (i  *  1 , . . . ,n) 
and  let  Y^  he  an  independent  normal  rv  with  EY^  -  0  and 
Var(Yn)  =  (-p)~l(l-p)(l+(n-l)o).  Then 


O', 


•V 


st 


100 


■Y„>IY0 


V- 


.+  Y  =0] 


Since  any  normal  density  is  PF^  it  follows*  from  Theorem  4.3,  that  T^ 
is  UK9  in  pairs. 

In  fact  we  can  obtain  a  stronger  result.  If  the  correlation  matrix 
of  T  is  of  the  form 


h 

n  \ 

/ 

(5.2.  1) 

• 

' 

- 

rn  J 

\ 

where  >  1, 

,  1  =  1 

,  .  .  . 

,  n  a  nd 

pa i rs . 

To  shot-/  it  note 

that 

every  it 

correlation  matrix  of 

a  multinomi 

1=1 


1/2' 


1/2 


-1 


((rrD1/2 . (rn-l)1/2) 


is  RR?  In 


be  a  sequence  of  independent  random  vectors 


■ft  x“\  e  = 

listributed  as  X_.  Clearly  j  is  a  multinomial  random 

i-1 

/er tor  with  correlation  matrix  (5.2.1).  Normalizing  Yv 


such  that  it 


•Ms  zero  means  and  unit  variances,  it  converges  in  distribution ,  hv  the 
mlti variate  central  limit  theorem,  to  a  multivariate  normal  random 
victor  with  correlation  matrix  (3.2.1).  By  Remark  (vii)  the  limit  in 
d i st r i hution  of  RR^  in  pairs  random  vectors  is  RR^  in  pairs.  The 
assertion  in  the  preceeding  paragraph  now  follows.  The  previous  result 
that  deals  with  the  symmetric  multivariate  normal  df  with  negative 
correlations  p  is  obtained  by  taking  r  *  1  -  p  in  (5.2.1). 

Every  multivariate  normal  random  vector  T  with  nonpositive 
correlations  is  NDS.  It  is  easy  to  verify  (2,4*)  directly  using  (1.2). 
Thus  T  satisfies  (2.3)  and  (2.3*). 


5.3.  Multivariate  hypergeometric. 

bet  )  have  the  probability  function 

i  n 


P(T  ■  . Tn  '  ‘„> 


M\  -1 
N/ 


M-  l  M  \ 

n  k 

1=1  1 

n  | 

n 

i«l 

iM"  l  h 

\  i-l  V 

>0,  l  t  O, 

i-1 


with  positive  integer-valued  parameter  vector  (N,  M_ , . . .  ,M  ,  M)  [see 

1  n 

.Johnson  and  Kotz  (1969]. 

The  multivariate  hypergeometric  df  is  the  conditional  df  of 
independent  binomial  rv’s  given  their  sum.  Thus,  by  Theorem  4.3  the 
hypergeometric  df  is  RR^  in  pairs  and  hence  it  is  also  EDS,  NUOD 
and  NLOD.  By  Remark  (iii)  the  loint  probability  function  of  (T, , . . . ,T  ) 
is  KR^  in  pairs.  A  special  case  of  this  fact  was  observed  by  Lehmann 
(1966),  p.  1144.  See  also  Ebrahimi  and  Ghosh  (1980).  By  the  discussion 
after  Theorem  4.3  it  follows  that  this  df  satisfies  the  RR^  condition 


of  Karlin  and  Rinott  (1980). 


We  will  show  now  that  T  is  NDS  by  proving  that*  for  0  <  t  <  N-l, 

—  —  n  - 


(r>.  3.  I  ) 


<TI . Vl>l  '„  ■  '„!  i  l<T . Vl>lTa  ■  *  11 


Then,  by  symmetry,  (2.4* )  follows.  Since  the  random  vectors  on  each  side 
of  (5.3.1)  have  multivariate  hypergeometric  df  we  actually  have  to  show 


0.3.2) 


U  %C  V 


where  U  *  (U^ , . . . #ll  ^ )  has  hypergeometric  df  with  parameters 

(N  +  1 ,  M  . . . ,M  ,  M  +  1)  and  V  »  (V-  ,  .  . .  ,V  . )  has  the  same  df  with 
1  n-I  —  1  n-1 

parameters  (N,  . ,  M)  where  N  *  N  -  t  -  1  and  M  *  M  -  t  -  1. 

1  n-l  n  n 

Denote  m  *  n  -  1 . 

Thinking  about  as  the  number  of  individuals  in  the  sample  in 

Category  i  (i  =  1 , . . . ,m)  it  is  easy  to  see,  by  conditioning  on  the  category 
of  t lie  first  individual  chosen,  that 

(u, . Um)  S-  <wJ*\...,W W1J£)  +  1,  . W<°)  if  the  first 

individual  is  in  category  £,  l  ■  l,2,...,m,  and  that  (U.  , .  .  •  fU  )  S"  (V_,...,V 

l  m  1  m 

if  tlie  first  individual  is  in  neither  of  categories  l,2,...,m,  where 
/ 0  )  (0  ) 

(W,  )  has  a  multivariate  hypergeometric  df  with  parameters 

1  m 

(N,  j,  M-l,  M  j , .  . .  ,M^,  M).  Thus  for  proving  (5.3.2)  it 

suffices  to  show  that 


0.3.3)  (w{t) . W^+l,  . W<°)  1  (Vx . Vm),  t  -  1 . m. 


We  will  prove  (5.3.3)  when  l  *  1 ;  the  proof  for  the  other  1 1 s  is 
similar.  Omit  the  superscript  1  and  consider  (W^,...,Wm)  * 
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which  has  a  hypergeomot  r  io  distribution  with  parameters 

<"•  :l.' . Vr  V  *>• 


Write 


M 


V-l 


P(v  -  v  ....  ,V  -  V  )  - 
11  mm 


M  ...M  M-  l  M 

i-1 


/ 


\ 


M  -  N 


[V1 . vm ’  N  j,Vi 

i  =  1 


/  V 


M1  -  vi » *  *  •  -  vm,  M  -  N  -  l  Ht  +  f  vt 

i«l  i=i 


/ 


Then ,  ft  is  easily  seen,  that  the  rv’s  , . . . tV  can  have  the  following 

interpretation:  First  divide  a  population  of  size  M  into  a  group  of  size 

N  and  another  group  of  size  M  -  N.  Next,  choose  at  random  individuals 

and  let  V  be  the  number  of  them  in  the  first  group  (of  size  N) ,  then 

choose  at  random  individuals  of  the  remaining  M  -  individuals  and 

let  be  the  number  of  them  in  the  first  group.  Continuing  this  way, 

finally  choose  M  individuals  of  the  remaining  M  -  M.  -  ...  -  M , 

J  m  i  m-l 

i nci ividuals  and  let  V  be  the  number  of  them  in  the  first  group. 

m 

The  rv's  W, , . . . ,W  may  have  a  similar  interpretation. 

1  m 

With  these  interpretations  it  is  easily  seen  that  if  the  first 
individual  chosen  in  the  *V  experiment.’  is  not  in  the  first  group  then 

V  W  V  (W.  +  1,  W„,...,W  ),  where  the  inequality  follows  from  (1.2). 
otherwise  (v.,...,Vm)  S-  (Wj  +  1,  w2»*“»wm>*  Thus*  unconditionally 

V  (Wj  +  l,...,Wm),  which  proves  (4.3.3)  when  9.  -  1. 
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3.4.  The  Dirich  let  distribution. 


bet  T  =  (T  . . . ,T  )  have  the  density 
I  n 


r(  l  V 

i=n  J 


i=o  J  "  en-1  n  V1  n 

f<ti . ‘n>  ■  v —  <*  - 1  v  ^ c,  •  cj  i  i  'j  i 1 

tt  r(o.)  j  i  j  j-i 


whore  the  parameter  vector  (0^,0^ * * *  *  * 0n>  satisfies  0j  ^1,  J  ■  0,1,. ..,n< 

The  Dirichlet  elf  is  the  conditional  df  of  independent  gamma  rv's 
given  their  sum.  Thus,  by  Theorem  4.3  the  Dirichlet  df  is  RR^  in  pairs 
and  hence  it  is  also  CDS,  NUOD  and  NLOD.  By  Remark  (iii)  f  is  RR^  in 
pairs.  Special  case  of  this  fact  is  Example  10  (iii)  of  Lehmann  (1966).  See 
also  Kbrahimi  and  Chosh  (1980).  By  the  discussion  after  Theorem  4.3  it 
follows  that  this  df  satisfies  the  RR^  condition  of  Karlin  and  Rinott 


(1980). 


To  prove  that  T  is  NDS  it  is  enough  to  show  that  for  t  <  t' 
“  n  —  n 


(5.4.1) 


<T1 . Vll|Ta  '  ?  1  <T1' 


then  (2.4*)  follows. 

It  is  easily  seen  that 


«-*n>-1|<Tl . Vl>lTo  ■  . Vl>lTn  ■  ^ 


Since  1  -  t  >  1  -  t*  (5.4.1)  follows  from  (1.4). 
n  —  n 


5.5.  Dirichlet  compound  multinomial  df. 

Let  T  -  (T1f..*,T  )  have  the  probability  function 
—  1  n 


3? 


»’(T)  =  t 


n:  r(  l  e  ) 

i  _n  J 


r(N  +  To.) 
J-o  1 


n  r(t1+fl1> 

1-1 


X 


t,  1  0,  l  t  <  N, 
J  j-1  J 


where  N  is  a  positive  integer  and  0^  >_  1,  j  *  0,1,..., n  [see  Johnson 
and  Kot z  (1969 ) 1 . 

The  Dirichlet  compound  multinomial  df  is  the  conditional  df  of 
independent  Pascal  (negative  binomial)  rv's  given  their  sum.  Thus,  bv 
Theorem  4.3  this  df  is  RR^  in  pairs  and  hence  it  is  also  CDS,  NUOD 
and  NLOD ,  and  by  the  remark  after  Theorem  4.3  it  satisfies  the  RR^ 
condition  of  Karlin  and  Rinott  (1980). 

To  show  that  T  is  NDS  it  suffices  to  prove  that 


'<T1 . Vl>lTn  ‘ 

then  (2.4*)  follows  by  symmetry.  But  the  random  vector  in  (5.5.1)  has 

the  Dirichlet  compound  multinomial  df  with  parameters  N  -  t^, 

f)  ,  ...,0  ,  .  Using  this  fact  it  is  easy  to  verify  (5.5.1). 

0  n-1 


5.6.  Df's  supported  by  negatively  tilted  surfaces. 

In  this  subsection  we  illustrate  two  simple  examples  of  such  df's 
ip\  The  general  idea  then  should  be  clear  and  will  not  be  detailed  here. 

het  (T^.T^T^)  have  the  uniform  df  on  the  surface 
f(t,,t2,t3)  1  (0,0,0):  tx  +  t2  +  t3  -  1).  Then  f (T, ,T2> | T,  -  has 


in 
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a  uniform  df  on  the  segment  —  (0>0):  t^  +  t2  =  1  “  C3^ 

which  is  clearly  stochastically  decreasing  in  t^.  By  symmetry  (2,4*) 
holds,  hence  (Ti*T2'T3^  is  NDS*  Clearly  ^Tj>T2,T3^  *s  a^so  CDS* 

Similarly,  if  (T^.T^T^)  have  the  uniform  df  on  the  surface  of 
the  unit  sphere  intersected  with  the  positive  (or  the  negative)  orthant 
then  (T1,T2,T3)  is  NDS  and  CDS. 

5.7.  NUOD  d f  *  s  in  reliability  theory. 

Buchanan  and  Singpurwalla  (1977)  considered  the  class  of  nonnegative 
multivariate  new  better  than  used  (NBU)  dfrs  that  satisfy,  for  all 

s  >  0,  t  >  0, 

(5.7.1)  PIT,  >  +  t, . Tn  >  ,n  +  tn)  <  P(T,  >  >  •„) 

'  P<T1  *  *1 . T„  >  *„>• 

It  is  well  known  that,  when  n  *  2,  (5.7.1)  implies  that  are 

negatively  quadrant  dependent,  that  is 


(5.7.2) 


P(T1  >  tl*  T2  >  t2)  iP(T]  >  tl)  P(T2  ’  t2)< 


Wo  will  show,  by  induction,  that  in  general  (5.7.1)  implies  that 

(T ...... T  )  is  NUOI).  From  (5.7.2)  we  know  it  is  true  when  n  -  2.  Assume 

I  n 

that 


(5.7.3) 


P(Tj 


•T 


n-1 


> 


t 

n- 


1 


) 


< 


n-1 

ir  P (T  >  t,). 
i-1  1  1 


Substitute 


to  obtain 
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